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We examine possible pairing mechanisms of superconductivity in PuCoGas based on spin- 
fluctuations or phonons as mediating bosons. We consider experimental data of specific heat C(T) 
and resistivity p(T) as input to determine a consistent scattering boson with the superconducting 
transition temperature of 18. 5K in PuCoGas . Irrespective to the type of boson, the characteristic bo- 
son frequency is found to be ~ 150K from the resistivity fitting. The spin fluctuation model is most 
consistent with the experimental resistivity, successfully explaining the anomalous temperature de- 
pendence ( ~ isoA'+T ) a ^ l° w temperatures as well as the saturation behavior at high temperatures. 
Assuming that the pairing state is non s-wave, the large residual resistivity pi mp ~ 20pQcm ~ 120 K 
suggests that an ideally pure sample of PuCoGas would have a maximum T c of 39 K. 

PACS numbers: 74.20,74.20-z,74.50 



I. INTRODUCTION 

Recently, superconductivity was found in PuCoGas at 
the amazingly high transition temperature (T c ) of 18.5 
K 1 . Considering the fact that the highest T c of f-electron 
based compounds 2 was ~ 2K , the T c of 18.5 K is an 
order of magnitude larger value than the previous high- 
est T c in f-clcctron based superconducting compounds. 
Therefore, the understanding the origin of this 18.5 K 
transition temperature in PuCoGas should not only pro- 
vide important information on the puzzling behavior of 
f-electrons but also shed light on the origin of the high 
transition temperature in cuprate superconductors. 

Let us briefly review the experimental data known 
about PuCoGas. First, from the specific heat jump at T c , 
the Sommerfeld coefficient is ^normal ~ 77 mJ/K mol. 
Second, from the T 3 phonon contribution in C(T), the 
typical value of the phonon frequency scale (0_d) is esti- 
mated to be 240 K. Third, the resistivity p(T) shows 
a typical S-shape behavior in its temperature depen- 
dence, which is often observed in spin fluctuating heavy 
fermion compounds such as UPt 3 [ 3 ]. Another impor- 
tant piece of information from p(T) is its magnitude; the 
value p(T = 300A') ~ 250^51 cm itself indicates strong 
scattering of the conduction electrons (an order of mag- 
nitude larger than the values of CeMIns (M=Co,Ir,Rh) 
superconductors having T c of ~ 2 K 2 ' 4 ). Forth, puz- 
zling is the data of the uniform susceptibility x(T)] it 
shows an almost exact Curie- Weiss temperature depen- 
dence of x(T) ~ 1/(T + T g ) with T e = 2K. This in- 
dicates that there are almost free local moments in the 
temperature range of 18.5 to 300 K with an effective lo- 
cal moment magnitude (0.68 ps) close to the value of 
the local moment of free Pu 3+ (0.84 ps) 5 - We believe 
that these local moments responsible for the observed 
x{T) are not coupled (or negligibly weakly coupled) to 
the conduction electrons and play no significant role in 



the transport properties as well as in the superconducting 
pairing. This does not mean that there are no interest- 
ing spin fluctuations - their contribution to x(T) might 
be much smaller than the local moment contribution or 
more probably their temperature dependence is not pro- 
nounced. If this is the case, the Curie-like x(T) should 
continue to exist below the transition temperature and 
should be observable if the diamagnetic part of x(T) is 
subtracted or suppressed below T c . 

As a possible pairing mechanism and pairing symme- 
try, we do not have much decisive data except estimates 
of a few energy scales. First of all, conventional phonon 
mediated pairing seems not unreasonable but only barely 
possible with values of Qd = 240 A, dimensionless cou- 
pling constant A = 0.5 ~ 1.0, and a typical value for 
the Coulomb pscudopotential p* = 0.1. The Allen- 
Dynes' T c formula 6 T c = g| cxp[- (^ffi+o^A)) ] pr0 " 
vides T C =16.7K for A = 1 and T C =2.9K for A = 0.5, 
respectively. On the other hand, although there is not 
yet direct experimental evidence, the existence of Pu 
f-orbitals participating in Fermi level crossing band(s) 7 
and the isostructure CeMIns (M=Co,Rh,Ir) compounds 
let us suspect the important role of spin-fluctuations to 
explain the normal state transport properties as well as 
the pairing mechanism in PuCoGas. 

In this paper, we examine two possible bosonic scatter- 
ing mechanisms, namely, phonons and spin- fluctuations, 
to consistently understand the available experimental 
data mentioned above. Our strategy is the following. 
Assuming each bosonic scattering, we try to fit the dc- 
resistivity data p(T) for its temperature dependence as 
well as its magnitude. From this fitting procedure, we 
extract the magnitude of the dimensionless coupling con- 
stant A and the typical energy scale of the corresponding 
boson. From these two numbers, we then can estimate 
T c using McMillan's formula. 
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II. FORMALISM 

We calculate the conductivity using the Kubo formula 
as follows. 

CT (T) = — V^fe) { —A 2 (k,u;)l 1 ] (1) 

K ' 3 ^ K } J AttT V ' ;L cosh 2 [ W /2T] J V ; 



where A(k, lo) = 2ImGn(k, to) is the one particle spectral 
density of quasiparticle in the conduction band and the 
retarded Green function of the quasiparticle is defined 
as Gji(k,uj) = 1 ,-» ; . All scattering- information 

is included in the self-energy S(fc, lo). Within the Born 
approximation, the self-energy is calculated as 



E(fc,w„) = <? 2 T^ / duo 



,B(q,uj') 
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iU n - lo' ito n + ifl„ - e k + q ' 



(2) 



where B(q, lo) is the spectral density of a bosonic propa- 
gator D(q,uo) = j 4j^r f^'^J, and g is the electron-boson 
coupling constant. After a Matsubara frequency 
summation, the imaginary part of lo + irj) is writ- 

ten as 



ImS R (T,uj + ir]) = g 2 N(0) 



du' 



(3) 



Q 

where n{u>) and f(ui) are the Boson and Fcrmion dis- 
tribution functions, respectively. N(0) is the density of 
states per spin at the Fermi level. Plugging the self- 
energy Eq.(2) into Eq.(l) and summing y\ , Eq.(l) gives 



He 2 

a(T) = —N(0) < v 2 > FS 



du> 



the overall magnitude. Finally, N(0) and < v 2 >fs are 
the values before they are renormalized by the bosonic 
scattering and, therefore, are difficult to be estimated 
from experiments. We rewrite Eq.(4) as follows. 

a(T) = ^ZN(0) < v 2 > FS (%[ 



4T L cosh 2 [w/2T] J /mS i? (T,< 
(4) 



A few remarks are in order for the above equations. 
First, the vertex correction is ignored. The justification 
is that when the scattering is local in space - technically 
meaning that the self-energy is momentum independent 
- the current vertex is not renormalized. This is consis- 
tent with the local approximation in calculating Eq.(2) 
and Eq.(3); consistent with this, we also neglect the mo- 
mentum dependence of the coupling g implying every 
quasiparticle is equally scattered by the assumed boson. 
Second, the self-energy is calculated only in the Born ap- 
proximation. Third, the Fermi surface (FS) anisotropy 
is neglected, resulting the factor | and the FS averaged 
Fermi velocity squared < v 2 >fs- With all these approx- 
imations, we should take the temperature power law of 
the calculated resistivity at low temperatures with reser- 
vation. Otherwise it induces an error of order 0(1) for 
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4T L cosh 2 [lu/2T] 1 ImZ R (T, lo) 
(5) 

In the above equation, iV(O) and v are the quanti- 
ties renormalized by the bosonic scattering and Z is 
the wave function renormalization parameter Z = 1 + 
dReS(u>) j duj. The above expression is very useful for our 
purpose. The renormalized quantities -/V(O) and v can 
be obtained from experiments. Furthermore, although 
Eq.(5) has an explicit dependence on Z, the implicit de- 
pendence of 7mS(T, to) on Z makes Eq.(5) a slowly vary- 
ing function of Z. The reason is because the real and 
imaginary parts of S(T, to) are related by the Kramer- 
Kronig relation, so that 9i?eS(T, to) /duo = Z(T) - 1 re- 
sulting M(T,w) - [Z(T) - 1]. 



III. RESULTS 
A. Spin fluctuations 

We choose the mean field type spin relaxational mode 

of B (i^) - [i ( T) +b tf%'r +m > > where r is thc ma §- 

netic Fermi energy 8 , I(T) = Iq + aT is thc parameter 
controlling thc distance from a magnetic quantum criti- 
cal point, and Q is a typical wave vector of the magnetic 
ordering. T ■ I(T) = co s f(T) defines the characteristic en- 
ergy scale of the fluctuations, b describes the dispersion 
of the collective mode of D(q, lo). The magnetic ordering 
vector Q can be in two dimensions or in three dimen- 
sions depending on the nature of the incipient magnetic 
-order. This dimensionality of magnetic ordering would 
) affect the power law of the resistivity at low tempera- 
tures. We assume 2-D in our calculations, reflecting the 
band calculations 7 . The overall magnitude C of B(q,u) 
is combined with g 2 N(0) in Eq.(3) to determine the over- 
all magnitude of £(T, to). This overall magnitude is de- 
termined once (Z(T) — 1) is fixed. Therefore, we do not 
need to determine separate values of g 2 , N(0), and C. 
Most importantly, Iq determines the low temperature 
T Q region of the resistivity before the inflection point; 
a ~ 4/3 was extracted from experiments 1 for T c < T < 
50K, but a dimensional counting for thc spin fluctuations 
model in 2-D gives p(T) ~ j~TaT a ^ ^ ow t em P era tures. 
The temperature variation of aT in I(T) also controls 
the high temperature saturation behavior of p(T), since 
increasing value of I(T) for larger temperatures reduces 
the scattering rate. 

Fig.(l) shows the typical results from Eq.(5) with Z 
values varying from 2 to 8. To be quantitative, the ex- 
perimental values of JV(0)(<— 7 = 77mJ/K 2 mol) and 
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v(= 8.8 x 10 6 cm/sec) 5 are used. Given these two ex- 
perimental input values, there is no free parameter to 
adjust the overall magnitude of the resistivity. The de- 
tailed temperature dependence of p(T) is controlled by 
I{T). r • Iq = oj s f = 150 K, a = 1 are chosen for illus- 
tration. As seen in Fig.(l), with increasing Z, the sensi- 
tivity of p(T) to Z becomes weaker, and for all values of 
Z the calculated p{T) is smaller than the experimental 
one by a factor of 3 to 4 (see Fig. 2). In order to fit the 
experimental p exp {T), we tune the value of v, which is 
a rather rough estimate from the critical magnetic fields 
(H C 2(0)) in superconducting state of PuCoGas, while the 
value of j exp = 77m J '/ 'K 2 mol is more reliable. We also 
chose to use Z = 4.6 which is the ratio of ~f exp / 'jband 7 ■ 
The result is shown in Fig. (2) in comparison with the ex- 
perimental Pexp(T). Input parameters are io s f = 150K, 
a=l, 7 = 77mJ/K 2 mol, and v exp = 4.78 x 10 6 cm/sec. 
The residual resistivity pimp = 15/if2 cm is added to the 
theoretical result (The better fitting of the low temper- 
ature part only would give Pi mp = 19/aQ. See the in- 
set of Fig. (2).) The overall fitting is satisfactory from 
low to high temperatures. In particular, the saturation 
behavior at high temperatures is well reproduced with 
the temperature dependent I(T). This is an expected be- 
havior since increasing temperature should shorten the 
magnetic correlation length such as £ _2 (T) ~ {Iq + aT). 
Some discrepancy between the theoretical result and the 
experimental one beyond 250 K (Fig. 2) would be due to 
the failure of the simple relation £~ 2 (T) ~ (/ + aT) at 
such high temperatures. To obtain better fitting in the 
low temperature region, one needs to allow a modifica- 
tion of the high temperature part. The result is shown 
in the inset. It indeed reproduces the experimental ob- 
servation p(T) - T 4 / 3 I 1 ] for T c < T < 50 K, while the 
correct theoretical form of p(T) at low temperatures is 

T 2 

~ 150K+T- 
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FIG. 2. Theoretical calculation with Z = 4.6 (open blue 
circles) and the experimental resistivity p(T) (red solid line). 
Input parameters are ui 3 f = 150K, a = 1, 7 = 77mJ/K 2 mol, 
and v exp = 4.78 x 10 6 cm/ sec and pi mp = 15^iSlcm is added. 
Inset: closeup view of low temperature region. For better fit- 
ting, v exp = 5.28 x 10 6 cm/sec and pimp = Wpflcm are used 
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FIG. 3. Theoretical calculations of resistivity p{T) with 
FIG. 1. Theoretical calculations of resistivity p(T) for vary- Einstein phonon with 60 = 240K for varying Z values {Z 
ing Z values (Z =2,4,6,8). Input parameters are u s f — 150K, =2,4,6,8). Input parameters are 7 = 77rnJ/K 2 mol, and 
a = 1 , 7 = 77mJ/K 2 mol, and v e x P = 8.8 x 10 6 cm/ sec. v exp = 8.8 x 10 6 cm/ sec. 
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FIG. 4. Theoretical calculations of resistivity p(T) for 
Z = 2 with two different Einstein phonon frequencies: 
6 D = 240K, 7 = 77mJ/K 2 mol, and v exp = 6.1 x 10 6 cm/sec 
(open green squares); 9d = 150K, and v e x P = 5.6 x 10 6 'era/ sec 
(open blue circles). The experimental resistivity p(T) 1 is red 
sold line. Inset: closeup view of low temperature region. 



B. Phonons 

Fig. 3 shows the resistivity calculated with an Einstein 
phonon B(q,uj) ~ 6(u> — 0d) with Qq = 240K, which 
is the value obtained from the specific heat C(T) above 
Tc 1 , for varying values of Z from 2 to 8. The input values 
are 7 = 77mJ/K 2 mol and v(= 8.8 x 10 6 cm/sec). As in 
the spin-fluctuations case (Fig.l), with increasing Z val- 
ues the sensitivity of p(T) to Z becomes weaker, and the 
overall magnitude of ptheor(T) is smaller than p exp (T) by 
the factor of 2 to 3 (compared the values at T=100K). In 
Fig. 4, we show the tuned theoretical results of ptheoriT) 
in comparison with the experimental p exp (T). As before, 
in order to increase the overall magnitude of ptheor(T), 
we tune the renormalized Fermi velocity v. 

In Fig. 4, we see that the theoretical resistivity 
Ptheor{T) with 9r> = 240A' (green open squares) has a 
higher power law region over a braoder range of low tem- 
peratures relative to the experimental data. For overall 
fitting of the data, we add a large impurity resistivity 
Pimp = 61//f2 cm, which makes the low temperature part 
in clear disagreement with experiment. To fit the low 
temperature region better, we need to reduce Qd- The 
blue open circles are the result with 60 = 15CW. This 
happens to be the same bosonic energy scale as used in 
the spin-fluctuation model fitting. It means that the re- 
sistivity data reveal a characteristic energy scale of the 
scattering boson to be ~ 150 K, irrespective of the ori- 
gin of the boson. An impurity resistivity p% mv — 36/xf2 
cm is added. The inset shows a close-up view of the low 
temperature region. There is a clear deviation between 
p ph (T) ~ exp [-e D /T] and p exp (T) ~ T 4 / 3 [ 9 ]. 



The phonon model has a critical defect at high tem- 
peratures with any reasonable parameters. While the 
experimental p exp {T) show a clear saturation behavior 
for T > 160 K, phonon scattering results in an ever in- 
creasing T-linear resistivity for temperatures beyond a 
fraction of 9d- This saturation behavior in resistivity is 
a long standing problem in A-15 compounds 11 , transi- 
tion metals, and some heavy fermion compounds such as 
UPt 3 [ 3 ]. While there is no general mechanism to explain 
this phenomena 12 , we can make the following remark. 

With the spin-fluctuation model, this saturation be- 
havior is naturally explained by the temperature depen- 
dence of the magnetic correlation length £ _2 (T) ~ 7(T) 
(see Fig. 2). On the other hand, phonon scattering needs 
to invoke a separate mechanism to explain the saturation 
behavior. Recently, Calandra and Gunnarsson 13 , assum- 
ing lattice vibrations couple with the electron hopping 
integral (HI), have shown that the saturation resistiv- 
ity arises at high temperatures due to a cancellation be- 
tween an increasing phonon population and an increasing 
electron kinetic energy. Without a specific model and 
numerical calculations, this saturation behavior can be 
described phenomcnologically in a two parallel resistor 
model ("shunting model") as 14 



P~ 1 (T)=pT d \ al (T)+p- 



(6) 



where p id eai(T) = p e - P h{T) + p lmp , and p e _ ph (T) is cal- 
culated by Eq.(5). p m ax is the maximum resistivity de- 
termined by the /-sum rule 13 , but here determined by 
empirical fitting. 
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FIG. 5. Shunting model fits to the experimental resistiv- 
ity (red sold line). The low temperature fit (open blue cir- 
cles) is with Vexp = 4.9 x 10 6 cm/sec and p ma x = 850/j.r? 
cm. The high temperature fit (open green squares) is with 
Vexp = 3.6 x 10 6 cm/ sec and p ma .x = 380^if2 cm. 

In Fig. 5, we show the best fit results with this model. 
For the electron-phonon resistivity p p h(T), we use Z=2 
and $d = 150K as in Fig. 4 and the overall magnitude is 
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tuned with v exp . As seen, the simple shunting model can- 
not fit the whole temperature region with any parameters 
in contrast to the successful cases of A-15 compounds 14 . 
The implication is either that the shunting model is not 
good enough to fit the fine details or that phonon scat- 
tering is simply not the correct model for the resistivity 
ofPuCoGa 5 . 



C. T c 

Before we estimate T c , we need to consider the effect of 
impurities in the sample. The effect of impurities on T c 
depends on the gap symmetry. For an s-wave gap - in the 
case of phonon pairing - T c is not reduced by nonmag- 
netic impurities 15 in the first approximation. However, 
for an unconventional non s-wave gap - in the case of spin 
fluctuation pairing - T c is reduced by any type of impuri- 
ties. As noticed in the previous section, the experimental 
resistivity has a large residual value of pi rnp ~ 20/iQcm, 
which should be a serious pair breaker for non s-wave 
pairing. 

Assuming pi mp = 20/iQcm in the spin fluctuation 
model and plugging Z=4.6 and v exp = 4.78 x 10 6 cm/sec 
in Eq.(5), we obtain ImY,i mp = Ti mp = 116K. This 
scattering rate is far larger than T c =18.5 K. From the 
theory of Abrikosov-Gor'kov 16 , the transition tempera- 
ture is reduced by ~ 7rr ^ mp . But the original Abrikosov- 
Gor'kov formula should be modified by the mass renor- 
malization factor Z as 7rFi ^p/ z , With Z=4.6 and A T c 
=T c o— T c ~ 20 K, we estimate the transition tempera- 
ture of an ideally pure sample to be T c0 ~ 39 K. There- 
fore, it appears that if the gap symmetry of PuCoGas is 
unconventional, as in most of heavy fermion supercon- 
ductors, PuCoGas is another true high temperature su- 
perconductor. The important question is then whether 
the parameters of the spin fluctuation model obtained 
in the previous section can produce a transition temper- 
ature of ^ 39 K. For a ballpark estimate, we use the 
Allen-Dynes formula 6 for T c , 

_ <lo> 1.04(1 + A) 

Tc = T2<r cxp[ — W ] ' (7) 

assuming /i* = for non s-wave pairing channel. 

The important parameter is < u >= 
[/ °° du a 2 B(u)]/[J^° duj a 2 B(uj)/uj\, which defines the 
characteristic energy scale of the pairing boson. The 
nominal characteristic energy scale of spin fluctuations 
is Wgf = 150 K, but the above Allen-Dynes definition of 
< oj > would give several times uj s f because of the 1/uj 
long tail of the spin fluctuation spectra B(u>). In reality, 
the effective cut off energy scale should be between these 
values, say, 3 ~ 4 oj s f S - Without solving the strong 
coupling T c equation with details of the band structure 
and full dynamics of the spin fluctuations, we cannot 
tell the precise value of this. With this reservation for 



< w >, the Allen-Dynes formula with A = Z — 1 = 3.6 
and < uj >^ u> s f = 150 K indeed produces T c ~ 35 K. 

For the phonon case, nonmagnetic potential scatterers 
do not affect T c for an s-wave symmetry gap 15 . There- 
fore, we do not need to consider impurity effects on T c . 
The Allen-Dynes formula with A = 1, and p = 0.1 gives 
T c = 10.44 and 16.7 K for < w >=150 K and 240 K, 
respectively. However, if we assume that the total mass 
renormalization (Z=4.6) is caused solely by phonon scat- 
tering, the effective coupling A ~ 3.6 should be used, and 
it would be strong enough to produce T c > 20 K with 

< uj >=150 K. 



D. Natural radiation damage 

PuCoGas should have radiation-induced self damage 
as a function of time. A preliminary measurement 1 indi- 
cates a decrease of T c at a rate of ~ 0.2 K/month, due 
to a radiation damage, which appears to be a quite slow 
suppression rate at a first look. Here, we estimate a the- 
oretical T c suppression with radiation damage and com- 
pare it with experiments. Because there is no study of 
radiation damage in PuCoGas itself, we use results from 
Pu metal (J-phase) 17 . Each Pu decays into U and an a 
particle. Both particles collide with nearby Pu nuclei cre- 
ating so-called Frcnkcl pairs consisting of a vacancy and 
a self-interstitial of Pu ions. In the case of <5-Pu, the total 
number of Frenkel displacements from one Pu is 0.1033 
per year. Given these data, the estimated number of 
displacements per Pu per month is ~ 0.86 %. Assum- 
ing the unitary limit of the impurity scattering strength, 
the impurity scattering rate is given by T imp = , 
and therefore T imp ~ 2 K. For s-wave pairing, this little 
potential scattering should have no effect on T c suppres- 
sion. However, it is not certain whether the displaced 
Pu ions would behave as potential scatterers or magnetic 
scatterers. For non s-wave unconventional pairing, T c de- 
creases even with potential scatterers. With a modified 
Abrikosov-Gor'kov formula, A T c =T c0 - T c = * v ™*l z , 
we obtain AT C ~ 0.3 K per month, assuming the mass 
renormalization factor Z — 5, consistent with the exper- 
iment. Therefore, we conclude that the T c degradation 
with natural radiation damage is consistent both with 
unconventional pairing and also with s-wave pairing, if 
the displaced Pu ions behave as magnetic impurities in 
the latter case. 



IV. CONCLUSION: 

We have found rather specific constraints on possible 
pairing bosons to be consistent with available experimen- 
tal data for PuCoGa 5 . We examined two possible bosons 
as a common source for the normal state resistivity and 
the superconducting pairing. 
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For a spin fluctuation model, we found that the charac- 
teristic spin fluctuation energy of 150 K is consistent with 
the resistivity data. Theoretical calculations of resistivity 
produce a satisfactory agreement with the experimental 
resistivity: both for the anomalous power law (~ T 4 / 3 ) at 
low temperatures and for the saturation behavior at high 
temperatures. The same spin fluctuations would pro- 
duce unconventional superconducting pairing, such as d- 
wave symmetry. In this case, the large residual resistivity 
Pimp ~ 20pflcm acts as a serious pair breaker. To survive 
with T c =18.5 K after this pair breaking effect, the orig- 
inal T c o should be ^ 39 K. The Allen-Dynes T c formula 
can support this high T c o with < lo >~ uj s f. If this is 
the case, we have observed another unconventional high 
temperature superconductor after the cuprate supercon- 
ductors. Perhaps this large impurity effect explains why 
UCoGa5 (pimp ~ 20p.fl cm) is not a superconductor 18 . 
Taken together, we think that spin fluctuations are the 
most likely source of superconducting pairing as well as 
the normal state resistivity in PuCoGas. 

For the phonon model, we found that the most consis- 
tent phonon frequency to fit the resistivity data should be 
also 150 K, having a discrepancy with the estimate from 
the specific heat measurement (#d = 240 K). However, 
theoretical calculations of the resistivity with phonon 
scattering have crucial defects to explain the experimen- 
tal resistivity First, it is impossible to produce p e xp{T) 
~ T 4 / 3 with phonon scattering at low temperatures. Sec- 
ond, the saturation behavior at high temperatures needs 
a separate explanation - this high temperature feature 
in resistivity has no direct relation with the supercon- 
ducting pairing mechanism, though. Putting aside these 
defects and pushing the possibility of the phonon pair- 
ing, the merit of phonon mediated s-wave pairing is its 
insensitivity to large impurity scattering observed in the 
sample. Using a typical strong coupling constant value 
A = 1, it is difficult to achieve T c of 20 K. But, as- 
suming that phonon scattering is the main source of the 
large mass renormalization of Z ~ 4.6, the Allen-Dynes 
T c formula with < u> >= 150 K can easily support T c 
of ~ 20 K. In this case we find some similarity with the 
phonon mediated A-15 compound superconductors, such 
as NbsSn. However, comparing with the spin fluctuation 
model, we think that the phonon model is unlikely in 
PuCoGas. Some key questions to be answered for the 
phonon model are: (1) the experimental p(T) ~ T 4 / 3 for 
T c < T < 50 K [ 4 ] is not easily reconciled with a phonon 
scattering mechanism; (2) the absence of superconduc- 
tivity in UCoGa5 also needs an explanation if phonons 
are the pairing boson. 



V. ACKNOWLEDGEMENTS 

We thank M.J. Graf, J. Sarrao, and N. Curro for dis- 
cussions. Work at Los Alamos was performed under 
the auspices of the US DOE. Y.B. was partially sup- 



ported by the Korean Science and Engineering Founda- 
tion (KOSEF) through the Center for Strongly Corre- 
lated Materials Research (CSCMR) (2003) and through 
the Grant No. 1999-2-114-005-5. 



1 J. L. Sarrao, L. A. Morales, J. D. Thompson, B. L. Scott, 
G. R. Stewart, F. Wastin, J. Rebizant, P. Boulet, E. Co- 
lineau, and G. H. Lander, Nature 420, 297 (2002); J. L. 
Sarrao, J. D. Thompson, N. O. Moreno, L. A. Morales, 
F. Wastin, J. Rebizant, P. Boulet, E. Colineau, and G H. 
Lander, J. Phys. Condens. Matters 15, s2279 (2002). 

2 C. Petrovic, P.G Pagliuso, M.F. Hundley, R. Movshovich, 
JL. Sarrao, J.D. Thompson, Z. Fisk, p. Monthoux, J. 
Phys.: Condens. Matter 13, L337 (2001); R. Movshovich, 
M. Jaime, J. D. Thompson, C. Petrovic, Z. Fisk, P. G 
Pagliuso, and J. L. Sarrao, Phys. Rev. Lett. 86, 5152 

(2001) . 

3 GR. Stewart, Z. Fisk, J.O. Willis, and J.L. Smith, Phys. 
Rev. Lett. 52, 679 (1984). 

4 H. Hegger, C. Petrovic, E. G Moshopoulou, M. F. Hund- 
ley, J. L. Sarrao, Z. Fisk, and J. D. Thompson , Phys. Rev. 
Lett. 84, 4986 (2000). 

J.D. Thompson, J.L. Sarrao, LA. Morales, F. Wastin, and 
P. Boulet, unpublished (2003). 

6 W.L. McMillan, Phys. Rev. 167, 331, 1968; P.B. Allen and 
R.C Dynes, Phys. Rev. B 12, 905, 1975. 

7 I. Opahle and P. M. Oppeneer, Phys. Rev. Lett. 90, 157001 
(2003); T. Maehira, T. Hotta, K. Ueda, and A. Hasegawa, 
Phys. Rev. Lett. 90, 207007 (2003). 

8 P. Monthoux, A.V. Balatsky, and D. Pines, Phys. Rev. 
Lett. 67, 3448 (1991). 

9 Instead of using Einstein phonon spectra, using a Debye 
phonon spectra it would produce p p h(T) ~ T 5 law at low 
temperatures 10 , again deviating from the experiment. The 
high temperature behavior is the same with either phonon 
spectra. 

N. W. Ashcroft and N. D. Mermin, Solid State Physics, 
Saunders College Publishing, p 526 (1976). 

1 Z. Fisk and GW. Webb, Phys. Rev. Lett. 36, 1084 (1976); 
P.B. Allen, W.E. Pickett, K.M. Ho, and ML. Cohen, Phys. 
Rev. Lett. 40, 1532 (1978). 

2 P.B. Allen, Physica B, 318, 24, (2002). 

3 M. Calandra and O. Gunnarsson, Phys. Rev. B 66, 205105 

(2002) . 

4 H. Wiesmann, M. Gurvitch, H. Lutz, A. Ghosh, B. 
Schwarz, M. Strongin, P. B. Allen, and J. W. Halley, Phys. 
Rev. Lett. 38, 782 (1977). 

5 P. W. Anderson, Phys. Rev. Lett. 3, 325 (1959). 

6 AA. Abrikosov and LP. Gor'kov, Sov. Phys. JETP, 12, 
1243 (1961) (J. Exptl. Theoret. Phys. (USSR) 39, 1781 
(I960)). 

7 W. G. Wolfer, Los Alamos Science No. 26, p 274 (2000). 

8 UCoGas 1 has p(T) ~ T 2 at low temperatures and 7 ~ 10 
mj/mol K 2 . Within the spin fluctuations model, consis- 
tently UCoGas has a much weaker coupling and a larger 



G 



Io > 250K. 



7 



